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Abstract 


The unified processing and research of multiple network models are implemented, and a new 
theoretical breakthrough is made, which sets up two new theorems on evaluating the exact electrical 
characteristics (potential and resistance) of the complex mxn resistor networks by the 
Recursion-Transform method with potential parameters (RT-V), applies to a variety of different 
types of lattice structure with arbitrary boundaries such as the nonregular mxn rectangular networks 
and the nonregular mxn cylindrical networks. Our research gives the analytical solutions of 
electrical characteristics of the complex networks (finite, semi-infinite and infinite), which has not 
been solved before. As applications of the theorems, a series of analytical solutions of potential and 
resistance of the complex resistor networks are discovered. In particular, three novel mathematical 
propositions are discovered when comparing the resistance in two resistor networks, and many 


interesting trigonometric identities are discovered as well. 
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I . Introduction 


Many complex scientific problems can be simulated by the resistor network model, such as 
many electrical and non-electrical problems in the field of physics, engineering and mathematics. 
The progress of circuit theory not only promotes the development of integrated circuit and 
electrification science but also promotes the interdisciplinary development of natural science. 
Resistor network models are so important that the issues of various disciplines can be studied by 
simulating resistor network, such as conduction in anisotropic disordered systems [1], percolation 
and conduction [2], Anisotropy in electrical conductivity [3], Nonlinear localized modes in 
two-dimensional electrical lattices [4], Electric circuit networks equivalent to chaotic quantum 
billiards [5], photonic crystal circuits [6], Manifesting the evolution of eigenstates from quantum 
billiards [7], dynamical signature of fractionalization [8], the processing of hexagonally sampled 
two dimensional signals [9], topological insulators and superconductors [10], topological properties 
of linear circuit lattices [11], topological spin excitations [12], three-dimensional printed meshes 
[13], topological insulator surface [14], fractional-order circuit network [15], the mean field theory 
[16, 17], finite-size corrections of the dimer model [18], lattice Green’s functions [19-22], resistance 
distance [23], and so on. In particular, two important equations of Poisson equation and Laplace 
equation [24, 25] can be simulated by resistor network model [26]. In addition, a real plane network 
of graphene exists in the real nature. 

It is well known that calculating the equivalent resistance between two arbitrary lattice sites in 
a resistor network is always an important but difficult problem since it requires not only the circuit 
theory but also the innovative algebra. For example, when the boundary of resistor network is 
arbitrary, it is usually very difficult to obtain the exact potential and resistance of the complex 
networks with arbitrary boundaries. In fact, the boundary is like a wall or trap, which affects the 
solution of the problem. Therefore, the reality requires us to create new theories to accurately 
calculate the electrical characteristics (voltage and resistance) of the complex circuit network. 

Let's review the research history of resistor networks. In 1845 Kirchhoff established the basic 
circuit theory (the node current law and the circuit voltage law). After 150 years, Cserti [27] 
calculated the two-point resistance of the infinite network by Green’s function technique, which is 
mainly focused on infinite lattices, and some applications of Green’s function technique were 
obtained in later literature [28, 29]. In 2004 Wu [30] formulated a different approach (call the 
Laplacian matrix method) and derived the explicit resistance in arbitrary finite and infinite lattices 
with normative boundary (such as free, periodic boundary etc.) in terms of the eigenvalues and 
eigenvectors of the Laplacian matrix, which relies on two matrices along two vertical directions. 


Later, the Laplacian matrix analysis has also been applied to impedance networks [31], after some 
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improvements, several new resistor network problems have been resolved [32-34]. However, the 
Laplacian approach cannot apply to the network with arbitrary boundary since it is impossible to 
give the explicit eigenvalues for the arbitrary matrix elements (associating arbitrary boundaries). 
But the boundary condition is important since it is real case occurring in real life. 

In 2011 Tan pioneered a new technique for studying complex resistor networks [35], which 
now is called Recursion-Transform (RT) theory of resistor networks [26]. The RT method depends 
on one matrix containing one directions, which is obviously different from the Laplacian method 
which depends on two matrices along two directions. With the development of the RT technique, 
many problems of non-regular network with zero resistor edges have been resolved [36-45]. In 
addition, the advantage of the RT method is that all resistance results are in a single summation 
differs from the Laplacian approach gave resistance results are in the form of a double summation. 
Recently, the RT method have been subdivided into two forms: one form is the matrix equation 
expressed by current parameters [38-44], which is simply called the RT-I method; another form is 
the matrix equation expressed by potential parameters [26, 45], which is simply called the RT-V 
method. Summarizing the previous applications of the RT (including RT-I and RT-V) method, it is 
not hard to see that the previous researches of the RT method all rely on the zero resistor boundary, 
such as the globe network[36, 44] belongs to cylindrical network with two zero resistor boundaries, 
the cobweb network[26, 40] belongs to cylindrical network with one zero resistor boundary, the fan 
network[37, 45] belongs to rectangular network with one zero resistor boundary, and the hammock 
network[34, 43] belongs to rectangular network with two zero resistor boundaries. Obviously, how 


to study the complex network without zero resistor boundary by the RT method is a question. 


r--- 
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Fig.l. An arbitrary mxn resistor network with two arbitrary boundary resistors, where n and m 
are the maximum coordinate values of (n, m). Bonds in the horizontal and vertical directions are 
resistors r and ro except for two arbitrary boundary resistors of rı and r2. 
This paper developed a new technique and improved the RT theory to allow us to study 


arbitrary resistor networks without relying on zero resistor boundary, which can derive the electrical 


properties (potential and resistance) of the arbitrary mxn complex networks with complex 
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boundaries. Here we build two new theorems lead to large problems to be resolved. Our study 
shows the universal RT method is very interesting and useful to solve the complex network. We 
focus on researching the electrical properties (potential and resistance) of Fig.1 and Fig.2 on two 
complex mxn resistor networks with two arbitrary boundaries by the advanced RT-V method, which 
have not been resolved before. It is worth emphasizing that the non-regular complex networks with 
two arbitrary boundaries are the multi-purpose network model because it can produce various 
geometrical structure as shown in Figs.5-10. Thus a large number of problems of resistor networks 


will be resolved by this paper. 
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Fig.2 A nonregular cylindrical mxn resistor network, where n and m are the maximum 
coordinate value of (n, m), with the unit resistances r and ro in the respective horizontal and 
vertical (loop) directions except for two arbitrary boundary resistors of rı and 72. 


From the above analysis, professor Wu [30] was the first to give several accurate equivalent 
resistance formulas for the regular resistor networks by the Laplacian matrix method, for the sake of 
comparative study, here we introduce two main results of resistor networks from Ref.[30]. 


Case-1. Consider Fig.1 with z =7⁄ =Ņ is a regular mxn rectangle network, where n and m 


are the maximum coordinate values of (n, m), resistors r and ro are bonded respectively in the 


horizontal and vertical directions, the resistance formula for Fig.1 is 


r r 
R a(d d,)= x -x+ — |y - y| + 
nxn (disd) maS) 1 2 m MA 
2 m n LC, j cosy, +5)9, —C,, cosy, + OP T 
(m+Dn+D AS r'(1—cos@,) + 7, (1-cos8,) 


where C,, = cos, +54, , O=infm+) , ġ=jī/(n+1), d(x,y) and d,(x,,y,) are 


arbitrary two nodes in the network. 
Case-2. Consider Fig.2 with 7 =7 =", is a cylindrical mxn resistor network, where n and m 
are the maximum coordinate values of (n, m), resistors r and ro are bonded respectively in the 


horizontal and vertical directions, the resistance formula for Fig.2 is 
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_ 2 
R,,,,(d,,d,) = Af = yp BERN |, —x,| 


1 m n Cs + C, —2C C cos 2 (y, = y)0, 


oJ X2] 


+ 
(m+)Dn+Dae r,'\(1—cos20,) + r'(1—cos¢,) 


(2) 


where C, , = cos, +545, 0 =in/(m+l), $= ja/(nt)). 


The above results were found for the first time by Wu. Later Refs.[32-34] improved the 
Laplacian matrix method to make it applicable to regular cobweb and hammock networks. However, 
the improved Wu method still cannot resolve the resistor network with arbitrary boundary, such as 
the networks with two arbitrary boundaries of Fig.1 and Fig.2. In addition, the equivalent resistance 


in Eqs.(1) and (2) are in the double summation not in a single sum. 
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Fig.3 The resistor sub-network with the potential parameters. 


II. RT-V theory and Poisson equation 


Consider two kinds of complex mxn resistor networks of Fig.1 and Fig.2, where n and m are 


the maximum coordinate values of (n, m). Assume 4 (0,0)is the origin of the rectangular 
coordinate system, and denoting nodes of the network by coordinate {x, y}. Assume the electric 
current J goes from the input d,(x,,y,) to the output d,(x,, y,). Denote the nodal potential of the 
sub-network is shown in Fig.3, and expressing the nodal potential at d(x,y) by U,,,,(x%,y)=V.”. 


We will study the complex resistor networks in four steps. 
The first step, setting up discrete Poisson equation based on the sub-network of Fig.3. By 


Kirchhoff law (X.V, =0) to set up the nodal potential equations along the vertical direction, we 


achieve a discrete static field equation (or call Poisson equation) for any network, 


(Ay + hA WV =r PS > (3) 
where h=r/r,, and I =J (0,,,-,,,) Contains the input and output conditions of the current, 
AV =V -V9 +V® and AV! =Ve"—2W+Vo denote second order discrete equation, and 
when x,#x, Eq.(3) reduces to the discrete Laplace equation (AZ +hA}V{ =0. For the 
arbitrary network together with the upper and lower boundary conditions, by Eq.(3) we are led to 


V, V, > Via a rl, > (4) 


k+l 


=A 


m+1 


where V, and J, are respectively (m+1)x1 column matrixes, and reads 


T 
YE VO a UT 6) 
I, = Ey f ro , 1? ge the al (6) 
and A,,,, is the matrix built along the vertical direction. For Fig.1 and Fig.2, the A,,, is 
2+h+bh —h 0 0 —bh 

—h 21+h) —h 0 0 
A mai = a ” d ? (7) 

0 0 -h 21+h) —h 

—bh 0 0 —h 2+h+bh 


where b= w/t , and r, is the resistor between (x,0) and (x,m) in Fig.2, when b=0 (n =œ), 


the A 


in belongs to Fig.l; when b=1 (n=rņ), the A,,, belongs to Fig.2. The purpose of 


mel 
introducing r, is to express two different resistor networks uniformly. 

The second step, consider the boundary conditions of the left and right edges in the network 
of Fig.1 and Fig.2. Applying Kirchhoff's law (3°17, 'V, =0) to each of the left and right boundaries, 
we obtain two matrix equations of boundary conditions, 

hV, =1A,,.-(2—hE]V, . (8) 


hV, =[A„a —(2-A)EI, , (9) 


m+1 
where h,=7,/n, h,=1/n, E is the (m+1)x(m+1) identity matrix, matrix A 


Eq.(7). 


Equations (4)-(9) are all the equations we need to compute the node potential. However, it is 


is given by 


m+1 


impossible for us to get the solution of the above equations directly. Thanks to the RT theory of Tan 
that gave the matrix transform method [26, 38-40] and we create the new technique here. In the 
following we are going to give the transformation technology based on RT-V theory. 

The third step, creating matrix transform. Firstly, we work out the eigenvalue t, of matrix 
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A which is given by solving determinantal equation of det|A —tE|=0 Gust b=0 and 


m+1? m+1 


b=1), yields 
t, =20.+h)—2hcos 6, (i=0,1,---,m) (10) 
where 0, =(1+b)iz/(m+1), and b=0 for Fig.1, b=1 for Fig.2. 


Next to transform Eqs.(4)-(9) by the following approaches 


PraA nma = diag{to, tt, Pr (11) 
X, =P, V, or V, = (Pas) X, (12) 
where X, = [ x i eee” i . Assuming P is the row vectors of matrix Pm+i, such as 
R= Costi a ah (13) 
Thus, we multiply Eq.(4) from the left-hand side by Pm+1, we get 
Xp = 6 Xp) — Xe TI (6, 5,5 ath ya (14) 


where Eqs.(11) and (12) are used. 


Similarly, applying P, to Eqs.(8) and (9) , we are led to 


m+1 
hX =(t, +h, -2)X0, (15) 
BX = +h=2)x", (16) 

The above Eqs.(10)-(16) are all essential equations for evaluating the node potential. 


The fourth step, solving the matrix equations (13)-(16). Selecting sve =F —x,)rJ as 
i=0 


the reference potential, by Eqs.(14)-(16), we obtain after some algebra and reduction the solution, 
a [Ate =, J, (17) 
2-b\ 2 
where b=0 for Fig.l, b=1 for Fig.2,and x, is defined in Eq.(26) below, and have 


OF _ poe. 
xE Éran Peab oy (751) (18) 
(t -2)G, 


where C,,, È, GP are, respectively, defined in Eqs.(19)-(25) below. 


ks? 

The RT-V theory. The above method of establishing recursive matrix equations with voltage 
parameters, implementing matrix transform and obtaining the solutions of matrix equations is called 
RT-V theory. The detailed content of the RT-V theory (Recursion- Transform theory with potential 
parameters) can be found by the above four steps in Eq.(3)-(18). Here we give a summary of the RT-V 
theory on the consists of four steps: The first step creates a main matrix equation of potential 
distributions along the Y axis (such as the build of Eqs.(3)-(7)); The second step builds the 
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constraint equations (including boundary conditions) of nodal potentials (such as setting up Eqs.(8) 
and (9)); The third step diagonalizes the matrix relation to reduce the equations from two dimension 
to one dimension(such as converting Eq.(4) to Eq.(14) et al.); The fourth step figures out the analytic 
solution of the equations (such as the results of Eqs.(17) and (18)), then makes the inverse matrix 
transform by Eq.(12) to derive the analytical solutions of nodal potential and gives the resistance 


formula by ohm's law. 


Ill. Two theorems of resistor networks 


A. Several definitions 
In order to facilitate and simplify the expression of the solutions of matrix equations, we define 


several variables of C,, and 4,4 for later uses, 


Cpi = cos(y, a, De, ? Cy = cos(y, =; y)0, ’ (19) 


2, =1+h—hcos6, + f(+h—heos6,)—1, 


A, =1+h—-hcosd,—.{(1+h—-hcos,)’ -1. (20) 


with 
: b=0 for Fig.1 
6, =(1+b)iz/(m+)), ; (21) 
b=1 for Fig.2 
And define variables FO, a, pÊ and GP for later uses by 
F? =) -49/0 -4), AFP = Fake (22) 
Qs = AFD +(h, -DAF h,= 1/1. (23) 
l pË sero. terse 
Es ee aie diets : (24) 
ne E a a 2 if x 2 x, 
GP = FR +(h +h, -2)F? +h, -Dh -DFS . (25) 


The above definitions are applicable to the entire article. All of these definitions are meant to 
illustrate the following two fundamental theorems, and we always assume that the electric current J 


goes from the input d,(x,,y,) tothe output d,(x,,y,) in our entire paper. 


B. Two fundamental theorems 
Theorem-1. Consider the arbitrary mxn resistor networks of Fig.l and Fig.2 whose 


maximum coordinate value is (n, m). Then the potential of node d(x,y) in the mxn resistor 


network can be written as 


m 


1 = 
eS ae (26) 
C.F ye Ss. 


Vax (% y) = 
where (Cpo Gr) = ae and ¢,, is defined in Eq.(13), ¢,, is the conjugate complex of 
k=0 
Ç and X © is the solution of the matrix equation (14) together with the boundary condition 
equations. Formula (26) is a general formula which is suitable for any resistor network model. 
In particular, when selecting Èv = 5% —x,)rJ as the reference potential, the potential of 
node d(x,y) inthe mxn resistor networks can be written as 


X-X, 


V n(x% y) = 
m+l1 


2 —b m es 
pi —— YX C 27 
m +] 3 x y.i ( ) 
where X=(x,+x,)/2, &, is the conjugate complex of ¢,, (there be ¢,,=¢,, if ¢,, is justa 


real number), and x, is a piecewise function 
SX, (28) 


and X® is given by (18) which is the solution of equations (13)-(16). 
Theorem-2. Consider the arbitrary mxn resistor networks of Fig.l and Fig.2 whose 
maximum coordinate value is (n, m). Then the resistance between any two nodes d,(x,,y,) and 


d,(x,,y,) in the network is given by 


1 2 eee _ XE, 
(Sy oa) i=0 J ` 


where X® is the solution of the matrix equation (14) together with the boundary condition 


Rixn (d,,d,) = (29) 


equations, Formula (29) is a general formula which is suitable for any resistor network model. 


In particular, for the networks of Fig.1 and Fig.2, the resistance between two nodes d,(x,, y,) 


and d,(x,,y,) can be written as 


(i) (Fr 
imal, 2b ge KOE RUE, 


R _ (d,,d,)= 
ma (did) m+1 m+lia J 


(30) 


where b=0 is the case of Fig.l, and b=1 is the case of Fig.2, and X® is given by (18) which 


is the solution of equations (13)-(16). 
The above two new theorems contain a wide variety of geometric structure of the network 


model, which can produce many new results of potential and resistance, and can create new 
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mathematical identity (see the section 6). In the following, we are going to prove the correctness of 


two theorems. 


C. Proof of theorems 
Consider the mxn resistor network with two arbitrary boundaries shown in Fig.1 and Fig.2, in 
the introduction, we have built the key Eqs.(4)-(9) by the RT-V theory, and converted the equations 
to Eqs.(14)-(16), and derived Eqs.(17) and (18). Now we will work out the exact eigenvalues of 
matrix A„„ in Eq.(7). Eq.(10) can be derived by solving equation det|A,,—fE|=0, and then we 


need to consider two cases below. 


One is For Fig.1, substituting Eq.(10) into (11) with b=0 in A 


1/42 v2 2 


cosy, cosv,A, = cosv,0, 


we get the eigenvectors 


m+1? 


m+ 


, (31) 


cos vð, COsv,A, = cosy,6, 
where v, =k+1/2, and @ =iz/(m+1). By careful calculation, the inverse matrix can be easily 
obtained 
P; =-2 [P]. G2) 


where [ ]/ denote matrix transpose. 


Thus, the term ¢,,; appeared in Eq.(13) can be specifically rewritten as 
Ea 6 gH 24.05 ys) (33) 
E= 6, =c0s(y + 5)9,, (i21) (34) 


Another is For fig.2, substituting Eq.(10) into (11) with b=1 in A the eigenvector is 


m+1? 


obtained after some algebra and derivation, 


1 1 1 vee 1 
| 1 expG 0) exp(i20) =- exp(imd,) (35) 
m+l | + : . : 4 
1 exp(i@,) exp(i20,) =- exp(imé,) 


where 0, = 2iz/(m+1) (i=0,1,2,---m) . According to strict calculations, the inverse matrix reads 
E 1 1 exp(—i 0) iii exp(—i (A9) : (36) 


1 exp(—imO, ) Bi exp(—1m6,, ) 


Thus, the term ¢,, appeared in Eq.(13) can be specifically rewritten as 
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a Ors OM sym (37) 


Spi =expliv9,), Č, =exp(-iy8,) (38) 
We find that Eq.(32) and Eq.(36) can be rewritten as a unified form below 
an ae pie aM 


4 1 C10 Ou ai Ciz 


P oe Fa ? (39) 

anh g e 
a a a a 
where Cac.) = ae ae „and ¢ , is the conjugate complex of ¢ ,. 
k=0 4 a 
Using Eq.(12), we have V, =(P)" X, , expanding this matrix equation, then we get 
a [ez 4 xe a) (40) 
(Sk i? A) i ° x 


Equation (40) agrees with the formula (26) that we need to verify. 


Further, we get (¢ aan z mil by comparing equation (39) with equations (32) and (36). 


And when selecting Ae =5(% —x,)rJ , we have Eq.(17). Substituting Eqs.(17), (33) and (37) 
i=0 


into Eq.(40), then Eq.(27) can be verified immediately. 
Next, we verify Eqs.(29) and (30), by ohm's law, we have 


R nxn (did )= = LIV, y) -V Ya) (41) 

Substituting Eq.(26) with x={x,,x,} and y= {y,, y,} into Eq.(41), we therefore obtain Eq.(29). 

Substituting Eq.(27) with x={x,,x,} and y={y,, Y2} into Eq.(41), we immediately obtain 
Eq.(30). Thus, two theorems are verified. 

In subsequent section we consider applications of theorems to arbitrary lattices. In all 

applications, we stipulate all parameters in Eqs.(18)-(39) apply to all resistor networks, and denote 


the resistors along the two principal directions by r and ro except for resistors on the left-right 


boundaries, and the input and output nodes of current are respectively at d,(x,,y,) and d,(x,,y,). 


IV. Electrical properties of complex rectangular network 


A. Nodal potential of complex rectangular network 


Consider the non-regular mxn resistor network shown in Fig.l, where the maximum 
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coordinate is (n, m), selecting X V® = ee —x,)rJ as the reference potential, the potential of any 
i=0 


node d(x,y) in the finite and sem-infinite networks can be written as 


U ma% y) fa Xx- x, h < VAE em = ee 
J “mtl m+ ps (L= cos 8)G° Cpi» (42) 
x — m Ab 1 Mor ee Cu 
asad) = 2 EC, (43) 
m m+ i +h- Theos Oy" = 


where @,=i/(m+1), and C,,, f°, GP are, respectively, defined in Eqs.(19)-(25). For 
Eq.(43), there be no, x, x, œ withfinite x,—-x. 
In particular, when x, =x, (means the input and output nodes of currents are at the same 


vertical axis), formulae (42) and (43) reduce to 


Usd) hh N C, i 7G dG © 


= : ; 44 
T m+1 £4 (cos 8G?" = 
Um e g Cy ODC g, (45) 
J m+15 [4 +h-hcosð}-1 | 
Proof of Eq.(42). For Fig.1, substituting Eq.(34) with y, ={y,, ya} into (18), we achieve 
(i) 

Cig 

O Pes Crs Boul, rJ, (0<k<n) (46) 


' (t, a 
Substituting Eq.(46) and (34) into (27) with b=0, we therefore achieve Eq.(42). 
For proving Eq.(43), when noo, x, x, —>œ with finite x,—x, it can be got a limit by 


using Eqs.(20)-(25) 


ge Re x 
li MYX t 2 k = 47 
Mga (t, 7-8 (47) 


x0 


So, substituting Eq.(47) into (42) with n —> œ, we therefore verified Eq.(43). 

Formula (42) is a meaningful result because the network of Fig.1 is very complex and has not 
been resolved before, which contains a lot of different network models since the different boundary 
resistors can produce different geometric structures. Here several special applications of formula 
(42) are given below. 

Application 1. When 7 =n =h, Fig.l degrades into a regular mxn rectangular network, the 


potential of a node d(x,y) in the network is 


pÈ (i) 
Uy x — x, ie XVI „C Yİ BLASTA 
J m+ -cosb ) FË 


=1 n+l 


(48) 


where B° reducesto BO =AF°AF® 


n-X, 


In particular, when x, = x, potential formula (48) reduces further to 


U_ Gy 6 GC, -C dCi, oe 
mxn = b Ya Ys p a6 
J m+1 >, (l -cos 0)F® Pins (49) 


n+l 


Application 2. When h, =0(7, =0), Fig.1 degrades into a fan network as shown in Fig.4(a), 
where r and 7, are the respective resistors along longitude (radius) and latitude (arc) directions, 
and the resistor element on the outer arc is r, (an arbitrary boundary resistor). The potential of a 


node d(x,y) inthe mxn fan network can be written as 


m — RÖ 


x — D C C. 
U(x, y) = X x, r+ 2r Iva Jı XVX) Hio (50) 
J m+1 m+1 ms AF +(h, -DAF® yi 


n-1 


where we redefine B°. =F°a! (if x<x,) and BO. =F°af, (if x>x,). 
Please note that a non-regular fan network (the outer arc resistor r, is arbitrary) is a scientific 
conundrum, which has not been solved before. Ref.[26] has researched just the regular fan network 


(the outer arc resistor is r, = 7), but our formula (50) with zn =y is different from the result in 


Ref.[26] because two results depends on the different matrices along the orthogonal direction. 


D 


A 


(a) (b) 


Fig.4. Two resistor network models. (a) is a Fan network with an arbitrary boundary 
resistor r2; (b) is an arbitrary hammock network. 


Application 3. When 7 =r, =0, Fig.l degrades into a hammock network as shown in 


Fig.4(b), the potential of a node d(x,y) inthe mxn hammock network can be written as 


— (i) (i) 
U(x, y) _ X—X, a 2r 2 DON oe = wa nic 
J — m+l m+l£ F® yis 


SD 


where we redefine BY, =F°F®, Gf x<x,) and BY, = FOF”, (if x= x). 
In particular, when d,(0,y,) and d,(n,y,) are respectively at the left and right poles, the 


potential of Eq.(51) reduces to 


U(x, y) _ n—2x 
J m+) (52) 


Application 4. Consider the input current J is at d,(0,y,) on the left edge, and the output 
current J is at d,(n,y,) on the right edge, the potential of a node d(x,y) in the non-regular 


mxn resistor network of Fig.1 is 


mrak C GC. ; 
U(x, y) _ n—-2x r+ pie 2,n—-x “yi Lx Yad (53) 


J 2m+l) m+1&  (1—cos6)G le 


where œ, is defined in Eq.(23). 


In particular, when 7,=0 and h =1 (7 =n), the potential (53) reduces to 


U(x, y) _ n—-2x r+ 2r — fia c. 
J Wm+ mG AFO O 


(54) 
Application 5. Consider a regular mxn rectangular network of Fig.1 with 7 =r, =ņ, when 
d (0, y,) is on the left edge, and d, (n, y,) is on the right edge, the potential of a node d(x,y) in 
the rectangular network is 
Uy) n-2x „4h SAO ZA ig: (55) 
J 2(m+1) m+1 a- cos6)F® vi 


i=1 n+l 


In particular, when y, = y,, the potential equation (55) reduces to 


m G) _ (i) C C 
U(x, y) _ n—2x h AF. AF’ Yİ | (56) 


Jo Xm+1) m+1S Fo 1—cos 6, 
Application 6. Consider d,(x,,0) is at the bottom edge, and d,(x,,m) is on the top edge, 
the potential of anode d(x,y) inthe non-regular mxn resistor network of Fig.1 is 


y— m D =| i (i) 
U(x, y) _x-x i, Poe Y Bess CC... Ss 
J m+1 m+1 a (1—cos 6,)G® i yi 


where ®© is defined in Eq.(24), and C,; = cos(5 6). 
In particular, when x, = x,, the potential equation (57) reduces to 


Uy na S-Di 


a 58 
J ne 1 i=1 d- cos 0, )G 0,0 ~ y,i ( ) 

when x,=x,=0 and 7 =r, =r, the potential equation (57) reduces to 
~m+t+i > ECC, (59) 


a- cos) FO ap aS 


n+l 


Application 7. Consider d (0,0), and d,(n,m) are on two diagonal nodes, the potential of 


14 


anode d(x,y) inthe non-regular mxn resistor network is 


Uyi n—-2x 1 Ariat- (Dna? 
J E 2(m+1) di 1 > i= cos 6,)G® vil y,i > (60) 


where œ, is defined in Eq.(23). 
In particular, when z =n, = ņ, the potential of Eq.(60) reduces to 


U(x, y) _ n—2x , < AF —(— DAFP 
J ~ m+" eee die cos 6) F 0S yi? 


=] n+l 


(61) 


Application 8. Assume Fig.l is a semi-infinite oœooxn network, and m— œ but n ,x and y 


are finite. Consider d (0, y,) is on the left edge, and d,(n,y,) is on the right edge, when 


r =r, =ù, the potential of anode d(x,y) in the semi-infinite œoxn rectangular network is 


AF, -C — AF, C, )C, 
U on (2y) _ "a" r (AFC, IO 46. (62) 
J m 40 (1-cos 6@)F 


n+l 


where C, =cos(y, +5), F, =(Ak-A)[(A-A) with 2=1+h-hcos6+J+h—-hcosé) -1. Please 
note that these definitions apply to all such issues as appear below. Eq.(62) can be derived by taking 
the limit of Eq.(55). 

Application 9. Assume Fig.l is a semi-infinite œoxn network, where m—>co butn, x and 
y are finite. When r =0 and h,=1(7,=n) =(0,y,) and d,=(n,y,), taking the limit 


m—> œ in Eq.(54), we achieve the potential in a semi-infinite coxn network 


U oxn (X, y) = 2F F, x 
ahar CeO (63) 


Application 10. Assume Fig.l is an infinite 0ox0o network, but x,—x and y,—y are 


finite, taking the limit n —>œ in Eq.(43), we have the potential in the infinite rectangular network 


AC APC 

U gah y) = -f 227) 0 5 (64) 
J m METIE 

Notice that Eqs.(62) and (63) belong to the case of a semi-infinite network, while Eq.(64) belongs 


to the case of an infinite network. 


B. Resistance of complex mxn rectangular network 


Consider an mxn rectangular network with two arbitrary boundaries shown in Fig.1, where the 


maximum coordinate is (n, m). Defining £? = 2°, , the resistance between two nodes d, (x, y,) 


and d,(x,,y,) in the finite and semi-infinite networks are respectively 


(i) 72 (i) (i) 72 
|x = x| FE KH &XFi Cpi = 2p CL Gi F PoC; 


m+1 mt+1i5 (1 -— cos 0,)G® >, (65) 


R „„(d„d,)= 


mxn 


2 2 7 ama] 
R (d,d,)= |x =| r+ Gs oe — 24, CO, Gi f (66) 
mxo "19 42 m+1 m+1i5 Jd+h-heos6) -1 


where @,=iz/(m+1), and C,,, GP are, respectively, defined in Eqs.(19)-(25). For Eq.(66), 
there be no, x, x, —>œ with finite x,—x,. Eq.(66) can be derived by taking the limit 
n—oco in Eq.(65). 

Proof of Eq.(65). For Fig.1, substituting Eq.(42) with k= x,,x, into (41), then Eq.(65) is 
verified. 

Eq.(65) is an exact expression which still contains a variety of resistance results with all kinds 
of boundary conditions because the left and right boundaries are the arbitrary resistors. For clearly 
understanding formula (65), we set h and h,, m orn as special values, and give several special 
cases to understand its application and meaning. 

Case 1. When h =1, the network of Fig.1 degrades into a rectangular mxn network with an 
arbitrary right boundary, then formula (65) reduces to 
n & Bury ~ 2BiCy yi + B2Co 


a roe 
r+ ; 
m+1  m+1& (1-cos0)[F® +(h-D)FO] 


n+l 


R madd) = 


mxn 


(67) 


where 8 reduces to {()) =AFO TAF. +(h, -DAF 2. 4]. 
Case 2. When h,=h,=1, the network of Fig.l degrades into a normal mxn rectangular 


network, then formula (65) reduces to 


(i) 2 (i) (i) 72 
|x > x| pa H 5 Ba Cpi T 2 Piz Cy Gi + PoC 
m+1  m+1 (d-cos 8) FË 


n+l 


Rma (dod,) = 


mxn 


l (68) 


where 2.) reduces to £') = AF\?AF,\”. . This problem has been researched in Ref.[30], and gave 


k,s 
Eq.(1) with a double sums. Clearly, our result (68) is different from Eq.(1). Two different results 
will be compared in the last section. This also shows that the equivalent resistance can be expressed 
in different forms. 


Case 3. When h =0, the network of Fig.1 degrades into a non-regular mxn fan network with 


an arbitrary boundary as shown in Fig.4(a), by Eq.(65), we obtain the resistance of a fan network 


D2 (i) (i) 72 
|x = x| 2r m-l Pi Ca = 2 piz C, Cya + ON 


m+1 m+f AF + (h, -DAF ® 


n-1 


R (d,,d,)= 


mxn 


, (69) 


where ;") is re-definedas Ë = FO[AF®, +(h,-DAF.?._,]. 


In particular, when h,=0, h,=1, the network of Fig.l degrades into a normal mxn fan 
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network, then formula (69) reduces to 


(i) (i) 72 
Pa —xX. m i —2 C ic tr CG. Í 
R (d,d )=! 2 i 2r aA X Piz Yı V Poa Y2» , 


70 
mxn m+1 ie i AF ( ) 


where 2?) is re-defined as £7’) =F, AF”? . This case has been researched in [37], but the result 


is different from Eq.(70), however they are equivalent to each other. The reason is that they choice 
the different matrix along different direction. This also shows that the equivalent resistance can be 
expressed in different forms. 


Case 4. When h =0, h,=2, the network degrades into a mxn fan network with double 


resistor edge, then formula (65) reduces to 


(i) (i) (i) 72 
poznal, 2r BBCi 2Bi2©y Cys + Bron m 
Ren Guede) = myi tmi Al + A," i (71) 


where ©) isre-defined as pË = FOU?“ +4"). 
Case 5. When h =h, =0, the network of Fig.1 degrades into a hammock network, so, formula 


(65) reduces to 
g Bapa 2B Cys + B2C3,s 


Rpa (dnd) = E24 2e — 


12 
m+1 ae Fo , Ue) 


Riza 


where B®? is re-defined as 8° = | ag ae . This problem has been researched in Ref.[34], but the 


result is different from Eq.(72)], the reason is that they choice the different matrix along different 
direction. 


In particular, when d,(0,y,) and d,(n,y,) are at the left and right poles, Eq.(72) reduces to 


Riven LO, y1}, {7, y=: (, = =0). (73) 


Case 6. When two nodes are on the same vertical axis, from (65) we have the resistance 


between two nodes d,(x,,y,) and d,(x,, Y2) 


(i) 2 
n 2 BO Cys En) 
Ryn {X> yib {x y) = m+l1 3 3) i cos 6, ` (74) 
Where 2 and G® are defined in (24) and (25). 
In particllar, when 7 =n = ņ , formula (74) reduces to 
AFM AF. (Cpa “Ca 

Rna (dX; yh {x y) = 5 F® ) 1- cos 0, (75) 

Case 7. When both d,(0,y,) and d,(0,y,) are at the left edge, formula (65) reduces to 

_ fi < AE +(h, - AFC (Cpe Cy, D 

R „x ({0, y,},{0, y} g G® 1=cos6, å (76) 


In particllar, when h,=1, d,=(0,0) and d,=(0,m), formula (76) reduces to 


n af, FO+h-DFo i 
R,„„({0,0},{0,m}) =— iž(1- FO 7 > Fup lcot’(} 8). (77) 
a n+l n 


It can be found that Eq.(77) agrees with the result in Ref.[41], clearly both of these results are 
verified indirectly by each other. 


Case 8. When both d (x,y) and d,(x,,y) are at the same horizontal axis, from (65), we 


have the resistance 


+1 m+l1 


i i i 2 lyg 
= |x, =z] H &¥ i -2893 + £53 cos [(y + z) | 78 
Raad yh {Xo y= m r+ l a 1—cos@ > (78) 


where £ and G® are defined in (24) and (25). 


Especially, when both d,(x,,0) and d,(x,,0) are at the bottom edge, Eq.(78) reduces to 


(i) (i) (i) 
xX, —X. r. m -2p £ 
E EA aE Ga 


r+ —— 5 TE 


ae 
mol mlg Jo (54): ©) 


When n =7 =Ņ, d,(0,0) and d,(n,0) are two corner points at the bottom edge, Eq.(79) 


reduces to a neat result, namely 


m( AF —1 
Ryu (O01 (m0 = ig rei SL a Jorgan: (80) 


n+l 


Case 9. When d,(0,y,) is on the left edge and d,(n,y,) is on the right edge, formula (65) 


reduces to 
Ra (0, 51), (2, yo} = BE 4 i E PCr PAC Cas Haars (81) 
m+l m4+15 (1—cos6,)G”” 
where a? =AF +(h, -DAF“. 
In particular, when h, =h, =1, Eq.(81) reduces to 
2 2 © 
Ry (Ory Toft yD = EG + Hh A, — Lori, (82) 


n+l 
Case 10. When d,(x,,0) is at the bottom edge and d,(x,,m) is on the top edge, then formula 
(65) reduces to 


n a pÈ -U-I BO + BY 


Rpa (0) rm) = ly 22 cot"(54)). (83) 


+ 1 m+ 1 i=l 2 GP 
where 2È and G‘ are defined in (24) and (25). 
Case 11. When d,(0,0) is at the coordinate origin, and d,(x,y) is an arbitrary point, by 


Eq.(65), the equivalent resistance is 


(i) 2 (i) (i) 2 
H m haC, i 2M, 0 nx), ¿Cni + 22C 


= KL Yast 
R nen 0, OF, {%, y) = mil mele (1—cos8,)G” ? (84) 


where ay,” and Ý are defined in Eqs.(24) and (25). 


In particular, when h,=h,=1 and y=m, Eq.(84) reduces to 


(i) i (i) (i) (i) 
—_ Xr h S AF, -2D AF, + AF. AF x 2 1 85 
Reon {0,0}, {,m)) = at 2F® cot (5 A), ( ) 


Case 12. When d,(0,0) and d,(n,m) are two diagonal nodes, by (65) we have the resistance 


between two maximally separated nodes 


hy m hag, = 2(-1)'h,h, + ha? 


—__Nr ln 2 1 
Rea OO m= 27 + OY TE co? (8), (86) 
where a? =AF +(h,-DAF®, 0. =in/(m+1). 
In particular, when 7 =n =, Eq.(86) reduces to 
ee! fy @) AB GN! A 
R a {0,0}, {2,m}) ehil Fa cot (54): (87) 


Please note that Eq.(87) is a desired equivalent resistance between two maximum separated nodes 
in an arbitrary mxn resistor network. This is an interesting result because it is simple and easy to 
research the asymptotic expansion for the maximum resistance. Refs.[46, 47] studied the asymptotic 
expansion by making use of the result (1). Obviously, the concise Eq.(87) is more conducive to the 
study of the asymptotic expression of the maximum resistance. In addition, there are similarities 
between equations (80) and (87), but only minor differences, where Eq.(80) is the resistance 
between two corner points at the bottom edge, Eq.(87) is the resistance between two maximum 
separated nodes on the diagonal line. 


Case 13. When n>, with x,=x,=0,and y,,y, are finite. Taking the limit of Eq.(76), 


we have the resistance in the semi-infinite network 


= n Š A, -1 (Cpi =C 
RelO 903) i AEG) 1- cos 6, $ (88) 


Obviously, case 13 is a semi-infinite network problem. The reason why the equivalent resistance is 
independent of the right boundary is that the network is infinite on the right. 
Case 14. When mn, h,=1,and y,,y, are finite. It means that the network is finite at 


the bottom and left but infinite at the top and right. Taking the limit of Eq.(88), we have 


r ez LCOS(y, + 39 —cos(y, + yor 


J, 1—cos0 sade’, (92) 


R x ({0, yı} {0, y} = 


0 
T 


where 2 =1+h—hcos@—,/(1+h—hcos 6)’ —1. This definition applies to all such issues as appear 
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below. 
In particllar, when m,n—o, and y,,y, >, h=1 (n=Ņ), but y,—y, is finite, taking 


the limit of Eq.(88), we have the equivalent resistance on the left edge 


al— — — 
Rese (0,94): 10,3) = 2 [E _ Fra. (90) 


Eqs.(89) and (90) belong to the problem of semi-infinite network, which has not been solved before. 


Case 15. When m,n, but x,—x, and y,—y, are finite, we have the resistance between 


two arbitrary nodes d,(x,,y,) and d,(x,,y,) 


or eae 7 |x| = 
Realy dy) = i| AO ag, (91) 
n zio [A+h-hcos0} -1 


Formula (91) can be proved by taking the limit m,n —>œ in Eq.(66). Eq.(91) shows that the 
equivalent resistance of an infinite network is independent of the boundary conditions. 

Notice that Eqs.(88)-(90) belong to the case of a semi-infinite network, while Eq.(91) belongs 
to the case of an infinite network. 

Since Eq.(65) is a profound and versatile formula that it contains multiple outcomes, for 
readers to better understand the meaning of Eqs.(65) and (90), we are going to again illustrate the 


meaning and usefulness by four simple examples below. 


r Bn 
ro r2 
Ao r A, r An 


Fig.5. A n-step resistor network with two arbitrary boundaries, the resistors in the horizontal 
and vertical directions are, respectively, r and ro except for two arbitrary boundaries. 


Case 16. When m=1, the resistor network of Fig.1 degrades into an 1xn resistor network 


as shown in Fig.5, the equivalent resistance between any two nodes in the resistor network can be 


written as 
r r. B, + 28 + L. 
R,,, (A, ,B )==—|, -x| + 2| “+ |, 92 
ica X a) 2 Ix, x| 4 G, ( ) 
r n| Ba 2b th 
R. (A,A yS— =x |) sl, 93 
me X i) 2 | 2 il 4 í G, l ( ) 
where £,,, G, are, respectively, defined in Eq.(24) and (25), and 
A=1t+h+ VW 4+2h, A=1+h—-Vh’ +2h. (94) 
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Fig.6. An arbitrary 2xn complex network of resistors with four arbitrary elements 
which are , respectively, 71, ro, r2 and r. 


Case 17. When m=2, the resistor network of Fig.1 degrades into an arbitrary 2xn resistor 


network as shown in Fig.6, where there be four arbitrary resistor elements (7, 7,7, 7), the equivalent 


resistance between any two nodes d,(x,,y,) and d,(x,,y,) can be written as 


h -269 + be i — 289 + (2) 
R,,, (A, + Ay) = A x,|+ a A Go 24 = 2) (95) 
o M4 gd e) 2) 4 B® 
r + 2B. f -2p 
Ry fA, B.) = 3 x, x,| + 2 aA GO 7 + 9GA = } (96) 
a) (2) (2) (2) 
r Bi +48, +4853 
Ry (Ay, > C, As sh gege (97) 


where 0<x <x, Sn, h, =r, /ņ(k=1,2 ), and È, GP are, respectively, defined in Eqs.(24) 


and (25). And F° =(4} -45 / (4-4), with 


=1Q+h+ V4), fy = 5 (24+ 3h + VOW 12h). (98) 


From the above derivation, we find that formula (65) is a generalized result, which is applicable to 
many network problems and summarized a variety of complex network models since it contains six 
arbitrary elements (}, r, f}, n, n,m). 

Case 18. Consider a regular mxn rectangular network, when m,n — œ (the network is finite at 


the left but infinite at the bottom, top and right.), by Eq.(90), the resistance on the left edge is 


T l+h 


1 
In particular, when h=1(r=n), Eq.(99) reduces to 


Rolt, y} (Oy +I) = 25. (100) 


Case 19. Consider a regular mxn rectangular network, when m,n — œ (the network is finite at 
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the left but infinite at the bottom, top and right.), by Eq.(90) we have the resistance on the left edge 


ReO yhy ad oo | h o pad it 
SO ares h+<(h 1) AS (101) 


H 


In particular, when A=1, from Eq.(101) we get 
({0, 0}, {0, 21) yn (102) 


Ros 


Eqs.(99)-(102) are two novel results which are given for the first time. 
V. Electrical properties of complex cylindrical network 


A. Nodal potential of complex cylindrical network 


Consider the non-regular mxn cylindrical network shown in Fig.2, where the maximum 


: ; ; a ae | ; Be 
coordinate value is (n, m), selecting XV re as the reference potential, defining 
i=0 


6, =2iz/(m+1), C,,-y = Cosy, —y)@,, the potential of any node d(x,y) in the finite and 


sem-infinite networks can be written as 


U (x, y) X-X h = P C y (Powe y 
mxn = try 1 i 2° 7 103 
J m+1 * Wma) > (1—cos0)G° a 
a P Ane _ zc 
U nmol% y) _xX—X, r+ r y i yy 1 22 (104) 


I mri Amda Jah hoso l 


where Ø“, and G; are, respectively, defined in Eqs.(22)-(25). For Eq.(96), there be n>, 


X» X, >œ with finite x, —x.Eq.(104) can be derived by taking the limit noo in Eq.(103). 
In particular, when x, = x, formula (103) and (104) reduce to 


U sra% y) = h S Cs — Cy,-y (i) 
J 2(m + 1) j=l (l —cos A)G° qvx? 


(105) 


m Cy -C 
Eo Li o (106) 
(m+D E /d+h—hcos6)— 


Proof of Eq.(103). For Fig.2, Ton P into Eq.(18), we achieve (O<k <n) 


XO = i 30 ~ rJ. (107) 
The substitution of (107) into Eq.(27) yields 
U(x, y) = x es h > n Br C yy ae 
J m+1 m+l1iz  21—-cos@)G° 
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rm & Be Sinl(y, — yel- BY, sin[(y, - yA] 
m+15 2X(1 — cos 8)G” ? 


4 (108) 


Because the elements 47, in the network is real number, the potential U(x, y) must be real number. 
Thus, extracting the real part of Eq.(108) to produce Eq.(103). 

Formula (103) is a meaningful result because the network of Fig.2 is very complex and has not 
been resolved before, contains a lot of resistor network models, where each of the different 
boundary resistor represents a different network structure. So Formula (103) can create many 
interesting results. In the following applications we always assume that the current J goes from 


d(x yı) to d,(x,,y,) except for special instructions. 
Application 1. Consider an arbitrary mxn cylindrical network of Fig.2 with 7 =r, =n, by 


(103) we have the nodal potential 


= © © 
U(x, y) 2X- h T ee Om — Le Oe 
J “m+ Xm+D4 (l-cos@)FO, ae 
where $°) reducesto 8°, =AF°AF®. . 
In particular, when x, =x, Eq.(109) reduces to 
UGy) n S Ga Gs po 
J Im +) HA cos VF Brave? (110) 


i=l 


n+l 


Fig.7. Two resistor network models. (a) is a cobweb network with an arbitrary boundary 
resistor r2, (b) is an arbitrary globe network. 


Application 2. Consider an mxn cylindrical network of Fig.2. When n =0, Fig.2 degrades 
into a cobweb network as shown in Fig.7(a), by (103) we have the nodal potential 


= @ © 
U(x, y) 2 XX, ack r uu a oe = ae Can 
J m+’ msl AFo+@,—DAF® ” 


(111) 


where 8° 


X,VX 


is redefined as B®. = Fa, (if x<x) and pÊ, =F°a) (Gf x2 x). 


In particular, when d,(0, y,) is at left edge, and d,(n,y,) is at right edge. Eq.(111) reduces 


to 
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Vey) pady pth g FP coso- 
J ~ 2m+1” m+ AF + (h, DAF?’ (112) 


Please note that the cobweb network with an arbitrary boundary has not been resolved before, the 
previous work only studied the normal cobweb network (the boundary resistor is r, =) [26], 
Eq.(112) is an original result. 

Application 3. Consider an arbitrary mxn globe network shown in Fig.7(b). That is to say 
that Fig.2 degrades into a globe network when 7, =r, =0, from (103) we have the nodal potential 


= @ @ 
U mL ¥- ph r $ Payal ee E oe 
J m+l1 ea = F® > 


(113) 


where we redefine BY. =FOF?, Gf x<x,) and BO. =F°F®, (if x2 x). 


n-X 


In particular, when d,(0,y,) is at left pole, and d,(n,y,) is at right pole, Eq.(113) reduces to 


UG, y) _ 2 
T mn vp 


Formula (114) is very simple and very interesting because the potential distribution is only related 
to the x and has nothing to do with y, which shows the nodal potential is equal in the same latitude. 
Application 4. Consider a non-regular mxn cylindrical network of Fig.2. Assume d,(0, y,) 


is on the left edge, and d,(n,y,) is on the right edge. by (103) we have the nodal potential 


U(x, y) h g Maan C,- y -haC = (115) 
J = tAmtDA d-cosG" 


where ap, = AF!) +(h,—IAF®, is defined in Eq.(23), and C, _, =cos(y, — y)6,. 
In particular, when y, = y,, Eq.(115) reduces to 
U(x, y) h m ha F hig (116) 
J = ary” n+) 5 (l—cos 8)G° 
when A =h, =1, Eq.(115) reduces to 
© © 
U(x, y)__n-2x r4 h s AF, Cy ï ZAF, C, Yr (117) 
J m+)" m+ FG (cos OF? ? 
when h,=h,=1, y, = y,, Eq.(115) reduces to 
m AFO _ Apo 
U(% y) _ ht AF, AF, G (118) 


J E7 : M + Xm +1) 2 (l-cos@)F°, ” 


Application 5. Consider a non-regular mxn cylindrical network of Fig.2. When the node 


d (0, y,)and d,(0,y,) are located on the left edge, Eq.(103) reduces to 
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U(x, y) _ ` =C CAA 
J -mie 12 ey, i etn 


where C, =cos(y, — y)@,, and a, is defined in Eq.(23). 
Application 6. Consider a non-regular mxn cylindrical network of Fig.2. When the node 
d(n,y,) and d,(n,y,) are located on the right edge, Eq.(103) reduces to 


U(x, y) _ wn Cy-y Cry Aa (120) 
J am 2 = o, GP’ 


One know the potential function have important application value for solving the Laplace 
equation. In this paper, the analytical solutions of node potential functions under various conditions 
are given, which provides a new theory for practical application. In particular, these simple 


equations of Eqs.(119) and (120) are very interesting and meaningful for applications. 


B. Resistance of complex mxn cylindrical network 
Consider a complex mxn cylindrical network with two arbitrary boundaries shown in Fig.2, 
where the maximum coordinate value is (n, m), and defining f° = pÈ. , The equivalent resistance 
between two nodes d,(x,,y,) and d,(x,,y,) in the finite and semi-infinite cylindrical networks 


are respectively 


id Alp -289 cos(y0,) + By 
Rsa (d»d,) = ET e 10088)" ; (121) 
m = 7 |x =x] 
Ce eles as iO, (122) 


m+1 0 m+lic' [+h—hcos6)—1- 
where @,=2iz/(m+1), y=y,—y,, È and GP are, respectively, defined in Eqs.(24)-(25). 
For Eq.(122), there be no, x, x, > with finite m. Eq.(122) can be derived by taking the 
limit n—>œ in Eq.(121). 

Proof of Eq.(121). For Fig.2, substituting Eq.(103) with k=x,,x, into Eq.(41), we therefore 
achieve (121). 

Formula (121) is an exact and exciting result because the network of Fig.2 is very complex and 
has not been resolved before, and contains a lot of resistor network models, where each of the 
different boundary resistor r, represents a different network structure. In particular, when taking 
some specific value for 7, and 1,, Eq.(121) gives rise to a series of special cases below. 


Case 1. Consider a non-regular mxn cylindrical network of Fig.2. When ņn=ņ, the 


resistance of Eq.(121) reduces to 
25 


ia =al h m 0 2B? cos(y8,) + 8È 
EF mn 4] © Ipro 
(m+1)& 1- cosO IF +(h, -)F°} 


n+l 


R xa (dod) = 


mxn 


l (123) 


where /;'’ reduces to pY) =AFP[AF®, +(h, -DAFP 1. 


Case 2. Consider a normal mxn cylindrical network of Fig.2 with 7 =n =ņ, the resistance 
of Eq.(121) reduces to 
k- TÈ in -2A cos(y0) + Br 
Tm +) (1-cos6)F® 


n+l 


R (dp d,) = : (124) 


mxn 


where ;’ reduces to fi? = AFOAF®. . 
Case 3. Consider a non-regular mxn cylindrical network of Fig.2, when h =0, the left 


boundary collapses to a pole, the network of Fig.2 degrades into a cobweb network with an arbitrary 


boundary resistor 72 as shown in far we have the equivalent resistance 


3 11 = 2B, cost yO.) + By3 
Eo AF® +(h,—1 DAF” 


n-l 


R (dy,d,) == 


mxn 


> (125) 


where /;') isre-definedas 8? = FO[AF®, +(h, -DAF À 


n-X, n-Xx, al i 
In particular, when h =0,h,=1, the network of Fig.7(a) degrades into a regular cobweb 
network, the resistance of Eq.(125) reduces to 


k- m RÒ _ 28 cos( y0) + pË 


1,1 
R,,,, (d,.d,) = ry va 


mxn 


, (126) 


where £;") is redefinedas gÊ = FOAF©. . 


Please note that case 3 has been researched in Ref.[37], but the result is different from 
Eq.(126), however they are equivalent to each other. The reason is that they choice the different 
matrix along different direction, where Ref.[37] set up matrix along the longitude, but this paper set 
up matrix along the latitude. 


Case 4. When h =h, =0, the left and right boundary collapse respectively to two poles, the 
network of Fig.2 degrades into an mxn globe network as shown in Fig.7(b), we have 


P- m RÒ Bi? cos(y8,) + By? 


R,,.,(disd>) = be jA oA, (127) 


mxn 


where £;') is re-definedas gË =FOF®,. 

Please note that case 4 has been researched in Ref.[36], but the result is different from 
Eq.(127), however they are equivalent to each other. The reason is that they choice the different 
matrix along different axes. This also shows that the equivalent resistance can be expressed in 


different forms. 
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Case 5. Consider a non-regular mxn cylindrical network of Fig.2, when d,(x,0) and 


d,(x,y) are on the same latitude, the resistance of Eq.(121) reduces to 


m (i) — 
R „a ({X,0},{%, y) = safi so), (128) 


m+liq GP 1- cos ð, 
where B° =a®a{_ and a is defined in Eq.(23). 
In particular, when h =h, =1, the network of Fig.2 degrades into a regular cylindrical network, 


the resistance of Eq.(128) reduces to 


i ARV ARS 1—cos(y6,) 
Ra ({x,0}, 1%, y}) = m+1 = Fo 1—cos@, (129) 
Especially, when d,(0,0) and d,(0, y) are on the left edge, Eq.(129) reduces to 
h 2(,_ f° \1-cos(y@) 
Re (O01(0.91) = 5S Fa) T= cos.’ i 


Case 6. Consider a non-regular mxn cylindrical network of Fig.2, when both d,(x,,0) and 
d,(x,,0) are on the same horizontal axis, we have 


(i) (i) (i) 
= x, =x n oh -28 +B 
R aren (x00) {%,,0}) = m+l aA (1 - cos 6,)G® ° (131) 


Case 7. When d,(0,0) is at the coordinate origin, and d,(x,y) is an arbitrary point, formula 
(121) reduces to 


(i) (i) (i) 
= ha,- 2ha n cosl yO.) + By 


2,n-x 


l ; 132 
i=l (1—cos 6, )G,” 


r 
Ren 10,0}, {5 y}) = ai 1 rt Gas i 


where œ}, and £,") are defined in Eqs.(24) and (25). 
In particular, when h,=h,=1 and y=m, Eq.(132) reduces to 


m AFO 2A + AF OAR 
i ous cos( yO ) x ni, (133) 


—_Xr r 
Ra OO} DoD = T t mt (1-cos8)G? 


Case 8. Consider a non-regular mxn cylindrical network of Fig.2, when d (0,0) is on the 
left edge and d,(n, y) is on the right edge, the resistance between two edges is 


m ha +ha® —2hh, cos( y8.) 
>t h han TMa I COS Y9) 
R a (10, 0}, {7, y} R aal 2(m+1) = (1-cos@,)G” ? (134) 


In particular, when h, =h, =1, Eq.(134) reduces to 


(135) 


m (i) 
__n h SAF. —00s(y4) 
Rma (10,03, {0 Y) = T t = (1—cos0,)F ý 


n+l 


27 


When h =0, Eq. (134) reduces to 


(136) 


i a a rh, m cg 
fem Gi ’ WMDs og! ts a4 = AF® + (h, -DAF d 


n-1 
Case 9. Consider a non-regular mxn cylindrical network of Fig.2, when d,(x,0) is at the 
bottom edge, d,(x,y) is an arbitrary node, and m,n—o, but x,y are finite, by(128) we have 


1—cos(y0)](1 + q2”) i 


z[ 
Rolt, Tao (137) 
where g=(A+h,-l)/(A+h,-l), Z=1+h-hcosO—J+h—hceos6y -1. 
In particular, when h =1, Eq.(137) reduces to 
Rolli Mae iy (138) 
When x—œ,but y is finite, Eq.(137) reduces to 
EU g (139) 


Ro({x,0},{x, y}) = a J+h=hcos6) -1 


Notice that Eqs.(137) and (138) belong to the case of a semi-infinite network, while Eq.(139) 
belongs to the case of an infinite network. 

From the above results we know formula (121) and (122) are two general results which, 
contains many results in a variety of lattice structures, can produce many new resistance formulae, 
and can create new identity (see Section 6). 

It is essential to take into account formula (121) again in order to help the reader further to 


understand its meaning, here two simple examples are given below. 


Fig.8 A 3D AxXn network with resistors r and v, in the respective horizontal 


and vertical directions except for rı and r2 on the left and right edges. 


Case 10. When m=2, Fig.2 degrades into a 3D Axn resistor network as shown in Fig.8. 


By Eq.(21) we have 6, =2iz/3 (i =1,2), substituting to Eq.(20) yields 


A =A, =1+3h+ a+ 3h)? =1. (140) 
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So we have ="), and 0,=27/3, 0,=42/3. By Eq.(121), we have the equivalent resistance 
between any two nodes 


-x P -289 cos(2ay/3) + By 
oy F Go 


k 
R,,(A, B) =! (141) 


where P, represents the nodes of A,,B,,C,, and 3"), G,” are, respectively, defined in Eqs.(24) 


and (25) . 
Eq.(141) is a general formula, which can produce two specific result below. 
Consider the resistance between two nodes A, and A,, there be y=0, Eq.(141) reduces to 


k-n]. 2n( 6-280 + 8% 
3/9 Go i 


(A, aye) (142) 


Raa 


Consider the resistance between two nodes A, and B, (or C,), there be y=1, Eq.(141) 


reduces to 


Ran 


k- D y ON (1) 
(A, ,B,)= Ry (A.C) = ae 1 aii 2,2 (143) 


Fig.8 is a simple and common network model, but getting the equivalent resistance has always been 
a difficult problem because of the complexity of the boundary conditions. Eq.(141) is given for the 


first time, which provides a new theoretical basis for practical application. 


Fig.9 A 3D OXn network with resistors r and y, in the respective horizontal and 


vertical directions except for rı and r2 on the left and right edges. 


Case 11. When m=3, Fig.2 degrades into a 3D oxn resistor network as shown in Fig.9. 


We can get the equivalent resistance between any two nodes by Eq.(121). By Eq.(21) we have 
0, =iz/2, substituting it to Eq.(20) yields 


A, =A, =1t+h+Jd+hP-1, 4,=14+2h+J0+2n)-1. (144) 


So we have £° =", and 6=2/2, 6,=7. By Eq.(121), we have 


he 250 needa £ Cy pa ies cos( yA, )+ Bo ies) 
16 ° 


R,,(Aps P, )= +2 GY GË 
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where P, represents the nodes of A,,B,,C,,D,, and p, GË are, respectively, defined in 
Eqs.(24) and (25) . 
Eq.(145) is a general formula, which can produce three specific result below. 


Consider the resistance between two nodes A, and A,, there be y=0, Eq.(145) reduces to 


-26% 1 0 (2) -2682 + (2) 
Rnb Ags A. )= a a gO : +2 6 GP = id (146) 
Consider the resistance between two nodes and B,, there be y=1, Eq.(145) reduces to 
k 
(1) 0 Ne (2) i 
ae; + +203 + 
Ratt) =hro4{ A GP 2) 3 G” 3 (147) 
Consider the resistance between two nodes A, and C,, there be y= 2, Eq.(145) reduces to 
(l) O, QQ) (2) (2) (2) 
k. nf Ba t22 + -23 + 
Rate bref > Gi = aT GO : (148) 


The case 11 tells us again that the general formula (121) is a meaningful and multipurpose result 


since just a 3D oxn resistor network has rich contents and many functions such as Eq.(145)-(148). 


VI. Comparation and Trigonometric Identities 


A. Proposition-1. A general trigonometric identity-1 


iT 
; b=. When m,n and x, Y, are 


ini = 1 __im 
Defining C; cos(y, +5)8, » and = a 


natural numbers, and O< x,,x,<n, O< y,,y, <m , we have the trigonometric identity 


ve Ic, į COS(x, + Dø, —C,, COS(X, + DA | 
E (1—cos8,) +h” (1-cosø,) 


2 


(i) (i) G) m2 

m C; =2 p3 CC t PC 

semti _yley lly, Yi ni 2,2 =“ yzi , (149) 
i=l (1—cos 6.) F) 


n+l 


where J) =AFOAF® and F? =A -79/4-7), AFP =F, -FP with 


A =1+h-hcos8, + Q+h-hcos8,} -1 , 


A, =1+h-hcos6.—/(1+h-hcos0,} —1. (150) 


Please note that the identity (149) is found for the first time by this paper. Identity (149) reduces a 
double sum to a single sum, which provides a new proposition and research method for 
mathematicians. 


Proof of the Proposition-1 
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Consider a regular mxn rectangular network shown in Fig.1 with 7=7, =ņ, Ref.[30] gave a 


resistance formula (1) by the Laplacian matrix method, which is in the form of double sum. 
However, this paper gives Eq.(68) by the RT-V method, where the condition and network structure 
agree with Ref.[30]. Obviously, the two results with different form in two different articles are 
necessarily equivalence because they are from the same network with the same coordinates. 
Comparing formula (68) with formula (1), we immediately obtain identity (149). 

We find Eq.(149) is an interesting identity for simplifying the double sum to be a single sum. 
In particular, when taking particular values of y,m,n and x,,x,, we have the following simple 


trigonometric identities. 


Deduction 1. When y, = Yy, = y, Eq.(149) reduces to 


1 1y4 | codia] 
y ng | costa, +4)4,—cos(x, +g, cos (y+ 5), 
n+l‘ 


22, (1—cos8)+h`'(1-cosø,) 


i i i 2 1 
2 B® -282 + BË | cos O +5)0, 


1,2 


@ 
al F 


n+l 
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1- cos 0, vey 


In particular, when x, =0, x, =n, Eq.(151) reduces to 


, 2 
m n -Cf cosd g) cosy +da, | m AF® 1 cos*(y +10, 
(1—cos6,)+h'(1—cos ¢,) > A 


n+l 


2 152 
n+l ae 1—cos 6, aa 


Deduction 2. When y, = y, =0, Eq.(149) reduces to 


| cost, + D4, — cos(x, + DA | 
n+l (l-cos8)+h'(1-cosø,) 


m (i) -2 p® + (i) 
= 1,1 1,2 22 2/1 
= 2 Fo cot“ (>0.). (153) 


n+l 


When m=l1, y=0, there be @ =7/2, Eq.(153) reduces to 


ig 155 7. | 
r | cost + Dh cos(x, + Dø, | _ nt = 263 + B33 , (154) 
mi2 1+h (1—cosg,) Fra 


J n+l 


where {2 =AFOAF® , FO =(Ak-2,5/A,-A) with A,=1+h+JhQ2+h). 
Deduction 3. When x,=x,=x, C, , =cos(y, + Da, , Eq.(149) reduces to 
2 2 1 
2 m n (Cpi =, a) COS («+ 5); 
n+1 i=l j=l ad —cos 6,) + h(i —cos¢@;) 
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m( AFOEAF® \(C, ,—C, 4) 
=—m+ nx | i hi 
=~ mt 1y, — T T -) en (155) 


n+l 


In particular, when x = 0, Eq.(155) reduces to 


m n (Cpi -Cy cos’ (Lg) +1 m Ee (C, = oe DH 
2 X 3 »2 i 274 = n(m My, y,| yh Yoi Yni , (156) 
n+l (1—cos,)+h (1—cos ¢,) n+l 


SEI 1 cose 
where AF\°AF® = F® — F and the following Eq.(176) is used. 


When y =0, y, =m, Eq.(155) reduces to 


pileos Iy T 
A we [1—( Df cost @,)cos(x+ dg, | 
n+l 


(1—cos6,)+h™'(1—cos ¢,) 


__mt1, CARD AFS, 
= ee Fo [1 (-1)'Joot”(5 6.). 


(157) 
n+l 
Deduction 4. When m=1, y, =0, y, =1, we have 8, = 7/2, Eq.(149) reduces to 
ie „ [cos(a +3), +0086 DAP BO +2 +p 2 (158) 
+15 o '(1-cos¢,) 2F n+l 
where F° =(Ak-2/(4,-2,) with 4,=1+h+Jh(2+h) and 7,=14+h-Ja2+h). 
In particular, when x,=x,=x and O< x <n, Eq.(158) reduces to 
7 1 i i 
2 æ COS (x+ 59; — AFOAFO | i (159) 
n+1jal+h(1—cos¢,) E n+1 


The above discoveries are interesting because they are found not in mathematics but in physics 


Obviously, according to the identity (149) we can derive a series of special trigonometric equalities 
when the different coordinates ( xi, yi) is made. 


B. Proposition-2. A general trigonometric identity-2 


Defining C, ; =cos(x, +50 , and 9, Si > $; -4 . When m,n and x,y, 


are 
natural numbers, and O< x,,x,<n, O< y <m, we have the trigonometric identity 

1 m 3 C j + Cc j =2C je J cos(2y0,) 

m+lizia (1-cos26,)+h`'(1-cosø,) 
2 (i) (i) (i) 
y n+1 2 fy 222 cos(2y0,) + 32 
i : 160 
Ha y) Zim +1) & (—cos20)F% = 


n+l 
where p =AFP AF! and FO=Ak-AN/A-A), AFO =F, -F with 
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2, =1+h—hcos26, +f(1+h—heos26,)—1, 


L 


J, =1+h—hcos26,—,{(1+h—heos20)° -—1. (161) 
Proof of the Proposition-2 


Consider a regular mxn cylindrical network shown in Fig.2 with ņ =r, =Ņ, we obtain a 


resistance formula (116) by the R7-V method. However Ref.[30] gave another resistance formula (2) 
by means of the Laplacian matrix method. Comparing formula (124) with Eq.(2), we immediately 
obtain identity (160). 
In particular, when setting special number values of y,m,n and x,,x,, we have the following 
identities. 
Deduction 1. When y=0, from (160), we have 
a tg Leost, +6, —cos(x, +DA P BO -2p + BO 


i l 22., (162) 
n+l (1-cos26,)+h'(1—cos¢,) — (1-cos26)F® 


n+l 


In particular, when m =1, we have cos20, =—1, then Eq.(162) reduces to 


ITIR Il 
ix | cosca + Ag, cos(a +D), | G BO -2p + pL 


n+1% 2+h '(1—cos¢,) 4F®? k oR) 
where Øo) =AFVAF)., FP =A -45/4-4), with 
A =142h+2Jhd+h), 7,=14+2h-2Jhn +h). (164) 
Deduction 2. When x, =x, = x , Eq. (160) reduces to 
> me cos’[(x+4)9,](1 - cos 2 y6,) 
m+ iia (1—cos20,)+h'(1—-cos¢,) 
` (= i> oes 2 aan Gave ) ne) 


Deduction 3. When m= y=1, we have cos20, =-1, 0, = ja/(n+1), from (160), we have 


fp et CAC, BO +2B2+ BY 2 (166) 
n+1ii2+h(1—cos¢,) 2F® n+l’ 


where F® =(4}-45/(-4), 4 and 4 are defined in Eq.(164). 


Please note that Eq.(166) is different from Eq.(158) because their 4, are different from each other. 


Deduction 4. When n=1, x=0, 6=iz/(m+1), 0< y <m, by (165) , we have 


“.1—cos(2y0,) _ B 
ieo = y(m+1- y). (167) 
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The above equations are given for the first time. 


C. Proposition-3. A general trigonometric identity-3 


Defining pË =AFOAF® , FË =t -49/4 -4), AA =1, 0 =iz/m, and 


n=X,? 


A, =1+h—-heos26, + /(1+h—hcos26,)’ -1, (168) 


and defining ¢=iz/(n+l), 444, =1, and 


u,=1+h'—h"'cos¢ + Jd+h'—h'cos¢y —1, (169) 


Assume natural numbers satisfy m21,n20, when natural number x, satisfies O< x, <m, 


and the arbitrary real number h>0O, we have 


n_ [cos(x +14 -cos(x +)gP 
1 X 2 ? 2 coth(mIn |, 
n+l JA+h”-h' cosg} -1 


m-1 -2p® 4 ie 
-h 1.2 

-x |+ , 
m” x| 2m <4 LSA (1—cos26, FO 


n+l 


(170) 


Proof of the Proposition-3 
Consider a normal m xn cylindrical network, where the maximum coordinate value is (n, m-1). 
When d,(x,,0) and d,(x,,0) are on the same longitude, by Eq.(131), we have 


x lpg t m- 1 Bi? - 289 i 9 
m ~ l= cos20)FË 


R({x,0},{x,0}) = bo (171) 


However, Ref.[48] gave another resistance formula (the parameters in Ref.[48] have been converted 


to be exactly the same as those in Fig.2) 


n x% | cost +D — cos(x, sba] 


— } > coth(m In 4), 
n+l J[O+h> -h cosg) -1 vu 


Thus formula (171) is equal to Eq.(172) since they have the same parameters in the same network. 


By Eq.(171) equals (172) to yield Eq.(170). That means, proposition-3 holds. 


R({x,,0},{x,,0}) = (172) 


Deduction 1. When x =0, x =n, ġ =iz/(n+1), @=ia/m, and u, is given by (169), 
Eq.(170) reduces to 
| xyi Lal 
n sin Gacostg | el (i) 
ape l cothtmin fig) =n Let (173) 
nA =A cA =l mae cos 20,) Fy 


where F? =(Ak-2,/(,-A), and A, is given by Eq.(168). 


Deduction 2. When m=2, ¢=iz/(n+1), 44 is given by (169), Eq.(170) reduces to 


i 
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l baT 
n [cos(x +5) — cos(x, Dg 0D 9B 4 
l 2 2 coth(In yi, ) = Ale, x,| + Bist 


(174) 
ntl Md+h'-hcos¢y -1 TA 


where pi) =AF”AF” , FP =Af-49/4-4). and A, is given by Eq.(164). 


Deduction 3. Since coth(ln y,) = (4; + 14)/( 4; — Z), substituting to (174) together with (169), 


which yields 
n 1 1 1 
h ap plea pin pafa E a 
where ø =izr/(n+1), {.) =AF”AF®, ‚and 4, is given by Eq.(164). 
Deduction 4. When m=1, ¢ =i/(n+1), Eq.(170) reduces to 
n [cos(x + ly —cos(x, + gp 
l 2 2 = |x, (176) 


n+14& 1—cos¢ 


where the following identity is used 


= -1 -1 2 
coth(In Ji= _ -m _ylth -h cosg) -1 (177) 


4+4 -2 h'(1-cos¢) 
We find that identity (170) is interesting because the left-hand side of the identity is the sum over n 
but the right-hand side is the sum over m, which provides a new mathematical identity for the 


application of mathematics. 


Vil. Conclusion and Comment 


This paper set up a universal RT-V theory (Recursion-Transform theory with potential 
parameters) and reveals the basic principle of electrical characteristics of complex resistor networks 
for the first time, such as two theorems of theorem-1 and theorem-2 are proposed, and the explicit 
electrical characteristics (potential and resistance) formulae of the complex networks are given, 
which contains the results of finite and infinite networks. 

It must be emphasized that previous theories (Mainly refers Green’s function technique and 
Laplacian matrix method) cannot solve resistor networks with complex boundaries, because the 
Green’s function technique is usually used to solve infinite network problems, and the Laplacian 
matrix method depends on the solution of two eigenvalues which relies on two matrices along two 
orthogonal directions. Using the RT-V method to study resistor networks just relies on one matrix 
along one vertical directions, which avoids the confusion of another matrix with arbitrary elements 
that cannot be solved explicitly, and also gives concise results in a single summation, such as the all 
equations given by this paper. 
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As applications of two theorems, the analytical solutions of the electrical characteristics 
(including potential function and equivalent resistance) in the complex mxn resistor networks with 
arbitrary boundaries are given, and many interesting results of the various types of resistor networks 
are produced. Please note that a non-regular mxn rectangular network (see Fig.1) contains arbitrary 
fan (see Fig.4(a)) and hammock (see Fig.4(b)) networks; a non-regular mxn cylindrical network 
(see Fig.2) contains arbitrary cobweb (see Fig.7(a)) and globe (see Fig.7(b)) networks. Thus the 
analytical solutions of the electrical characteristics given by this paper have general significance, 
which applies to a wide variety of lattice structures, and means the Laplace equation with complex 
boundary conditions is resolved by modelling resistor networks. The trigonometric identities given 
show that our research work established new research ideas and approaches for the study of 
mathematical identities. 

In addition, resistance formulae (65), (66), (121) and (122) et al. can be extended to impedance 


networks since the grid elements y, can be either resistors or impedances in Fig.l and Fig.2, For 


example, assume 
r=Z,=R+jolL, ņ=Z,=jl@C, (178) 
that we can therefore study the arbitrary mxn RLC network if we do a plural analysis [31, 41] to the 
resistance results obtained in this paper. 
Our resistance formulae can also apply to fractional circuits, for example, the frequency 


domain impedance of fractional capacitance and inductance are respectively [15] 


Z, = a Leos(B A + jo Lsin(B > (179) 
Za = l cos(a Z) — j ! sin(a Z) (180) 
© OC o to 8 


When we replace r and ro with Z, and Z,, and use the plural analysis [31, 41], the equivalent 


impedances of the fractional mxn RLC network with arbitrary boundaries can be obtained. 
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